Triplon Supersolid of Cold Atoms in a Ladder-Shaped Optical Lattice 
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We investigate the ground state phase diagram of the cold atoms in a ladder-shaped optical 
lattice by using stochastic series expansion quantum Monte Carlo method. We demonstrate how 
triplon supersolid could emerge, in which the triplon superfluid can coexist with the long-range 
triplon solid. In the triplon superfluid the quasi-condensation appears in the form of power-law 
decaying phase correlations. We suggest that a quasi-supersolid phase could be charted in the phase 
diagram, in which both superfluid and solid orders are quasi-long range ones. We also discuss how 
this phenomenon can be realized and detected in real experiment. 
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Introduction. — With the advantage of cleanness and 
controllability, the cold atoms system in optical lattice 
has been a powerful platform to simulate a variety of 
quantum system. Many works in experimental and the- 
oretical aspects have been made in the past few decades 
Various types of optical lattices with interesting 
geometric structure, such as triangular, kagome and hon- 
eycomb 0-@|: could be realized by adjusting the poten- 
tial formed by lasers. Further more, the interactions be- 
tween atoms can also be tuned by Fcshbach resonance 
fioj . With the aid of ultracold gas in optical lattice, the 
quantum phase transition from superfluid to Mott insu- 
lator is also observed [l|. Beyond that, we ask whether 
some fascinating quantum states, such as supersolid (SS) 
which is not observed in real materials so far, could be 
implemented by using cold atoms in some new kinds of 
optical lattices, such as a ladder-shaped optical lattice. 

Featured by the coexistence of both diagonal and off- 
diagonal long-range order, SS is an intriguing exotic 
phase of matter, of which research works date back to 
fifty year ago. Aroused mainly by the tentative exper- 
imental evidence for the presence of superfluid (SF) in 
solid 4 He although it is confirmed that all of those 
early results were misunderstood recently [l2j ]. an up- 
surge in quest for supersolid has been lasting for nearly 
a decade 13l42l|. In another route, lattice SS have 



been paid much attention, especially after the advances 



in optical lattice technology [lfj, l22|. Various schemes 
for lattice systems from one to three dimensions have 
been designed 23-3(|. Recently, the cold atoms in iso- 
lated double- well have been realized in experiment [H-Q . 
Could we go beyond those experiments to design a new 
system and realize a novel SS phase in optical lattices? 

In this Letter, we investigate SS phase by extending 
the isolated double- well |2|-l5[ to a spin dependent ladder- 
shaped optical lattice, as illustrated in Fig. [TJ It's can be 
viewed as double- well lattice that are coupled along y di- 
rection. Comparing with real compounds, the advantage 



of this system are the cleanness and highly controllabil- 
ity of the interactions which allow us extend previous 
studies of ladder system 31 - 3^1 to the parameters region 
where the triplon SS can emerge. This novel phase can 
be observed in ladder-shaped optical lattice by varying 
the lattice depth, the zeeman field and the interaction 
strength via the Feshbach resonance. When the interac- 
tions between cold atoms is larger than its kinetic enough 
i.e., U^>t, the system enter Mott insulator region. The 
dynamics of atoms just involves spin freedom of degree. 
By adjusting the zeeman field and the hopping of cold 
atoms with different spin, the triplon SF is induced by the 
power-law decaying phase correlations, while the triplon 
solid emerges as the result of triplon interactions. In a 
proper parameters region, triplon SF and solid appear 
simultaneously to form the triplon SS. 

The cold atoms in a ladder-shaped optical lattice. — A 
spin dependent ladder-shaped optical potential is con- 
structed as Fig. [T] We assume that only two of the in- 
ternal ground states participate in the dynamics, which 
we call spin up and down, | and | ],). By second-order 
perturbation theory for coupling of the ground state lev- 
els to the excited atomic states, the potential can be 
obtained as, U a = V x \, a cos 2 (k x x) + V X 2, a cos 2 (2k x x) + 



Vy.<j cos 2 (k y y), where ct=|, 4- represents the two states of 
cold atoms and V x \ j(T and V x i,a are the barrier heights of 
the two standing waves along the x direction, V yj(7 is the 
one along the y direction. In experiment a harmonic po- 
tential with frequency u> z should be added in z direction 
which restricts the moving of atoms along it. To simu- 
late a decoupled two-leg ladder potential, 4142, a- > V x \^ a 
is required and the ratio of potential heights between 
intra-ladder to inter-ladder along x direction should sat- 
isfies AV = (W x2 ^-V xha ) 2 /(AV x2t<T + V xha ) 2 < 1. The 
widths of the two barriers along rung direction and rail 
direction are a = a/k x and b = ir/k y , respectively, with 
a = arccos(I4i, (7/4142,(7)- With the considerations, then 
by using harmonic approximation in each bottom of trap, 




FIG. 1: (Color online) A ladder-shaped optical potential is 
formed by three standing wave lasers indicated as a pair of 
counter yellow arrows V x i,V X 2,V v . Two V x i,V X 2 of them 
along x direction and another one V y along y direction, where 
V X 2 (short yellow arrows) has twice the period of V x i (long yel- 
low arrows) . The cold atoms with two internal states denoted 
by spin up 1 1) (red up arrows) and spin down 1 4-) (blue down 
arrows) are trapped by this potential. An antiferromagnetic 
(AF) configuration with staggered up-and-down moments on 
the ladder is showed. The white lines in the potential contour 
lines represent the ladder. The effective strength of interac- 
tions is J± along x direction, Jjf, J| along y direction, h is 
zeeman field. 

this system can be described by Hubbard model as, 

h = ^2 (- W c i,<7 c j> + h - c ) + u n ^2 n *,t n a 

i,j,a i 
+ ^2,U a n i ^(n i ^ - 1) - h^2(rii^ - n a ), (1) 

i,a i 

where a zeeman term h is also introduced. c\ a (ci ir7 ) is 
the creation (annihilation) operator of boson or fermion 
in site i with spin a. The coefficients t^ tCr , with fj, = _L, ||, 
represent intraladder hopping along rung, rail direction, 
respectively. JTj-j., U a are the on-site interaction between 
atoms. We can work out those interaction parameters 
as a function of V X \ >(T , V x i,a: Vy,a-, a s , E xr , E yr , where a s 
is the scattering lengths and E x , yr = h 2 k 2 y /2m are the 
atomic recoil energy along x, y direction respectively. 

Considering t^ <C C/^, U a and 1/2 filling, by a second 
order perturbation theory, this system will be described 
by an effective Hamiltonian with spin freedom of degree 
only, H = ZiMSfS? ± J*(SfS* + - Sf, 

where S = ^c' ac, a is pauli matrix and = (cj, cj). The 
upper sign before is for fermionic atoms and down sign 
for bosonic one. The interaction coefficients for bosons 

2(i 2 -t-t 2 ) £ 2 i 2 

can be figured out as, J* — — J r- T — — -Jf- — -rr- , J?, = 

4 V t/4 ) h = 2h + ^ - For fermion, we need to 

omit the last two term in J* and h. 

Since the lattice is bipartite, without loss of generality, 
we focus on the region with positive of the spin XXZ 

model H and hereafter we take Jf = = J± = 1 as 
energy unit. The QMC simulation [341 ] is performed on a 
ladder with length L and temperature T = 6/L. A mean- 
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FIG. 2: (Color online) Ground state phase diagram with Jjf = 
0.35. The doted lines (Gray) are results from the mean field 
theory. The scatter lines (Black) are QMC results. There are 
four orders which are antiferromagnetic (AF), full polarized 
(FP), and triplon solid (TS) orders as shown in the insets. The 
triplon superfluid (SF) phase occupies a considerable area. 
The triplon supersolid (SS) phase is sandwiched between the 
AF and TS phase. The quasi-SS (q-SS) phase, which exhibits 
algebraically decaying diagonal TS order, is also charted. 

field (MF) theory is also employed for comparison. The 
ground-state phase diagram is illustrated in Fig. [5] 

The triplon superfluidity. — Let's start from the 
isotropic case, i.e., J* = J| = J||, which has been exten- 
sively studied in the context of real materials [35j ]. The 
four eigcnstates of jth decoupled rung provide very help- 
ful basis: the singlet state, \ s )j = ^(\ti)j ~ |lt) 3 -)' an d 
three triplet states, = Itt),-, l^^^dt^ + l^t),-), 
\t-)j = |44-)j- If 110 ( or weak) magnetic field is present, 
the ladder is in a gapped SD phase, whose ground state is 
roughly a chain of singlets. While the system will be in a 
FP phase in the strong magnetic field limit. The SD and 
FP phases can be readily expressed by MF wave func- 
tions, \iPsd)=Y\j \ s )j and IV'FpHIIf \ t +)j> respectively. 
If moderate field is applied, a gapless triplon SF phase 
emerges, which is induced by the quasi-condensation of 
tripon when the energy level of |£+) is down to the vicin- 
ity of state \s). And the upper \to) and |t_) bands 
could be integrated out from the renormalization point 
of view. One can introduce triplon bosonic operators 
b] = - fi+)that fulfill b] \s) j = \t + ) r Then the 

system could be mapped to a ID hardcore bosc-Hubbard 
model. It can be exactly solved by Bethe Ansatz or ap- 
proximately by Bosonization method, gives the Luttingcr 
liquid (LL) nature (3||. In this regime, the state |t+) un- 
dergoes a quasi-condensation characterized by a power- 
law decaying phase correlations (bjbj +r ) ~ r^ 1 / 2K cos 7rr, 
where K is the LL parameter. Although true condensa- 
tion is absent in ID systems, the triplon superfluidity is 
still possible due to the power-law decaying phase cor- 
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FIG. 3: (Color online) (a) The desity peak in momentum 
space (n w ), (b) The triplon structure factor S (n) for different 
system sizes, (c) The extracted index a and j3 versus ft for 
the isotropic case Jjf = J| = 0.35. (d) The staggered mag- 
netization m 3 and (e) its corresponding staggered magnetic 
susceptibility \s for Ju =0.35 and ft = 0.5. 
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FIG. 4: (Color online) (a) The density peak (n T ) and (b) 
the structure factor S(n) versus L in log-log scale. The 
parameters are: SF (ft = 0.96, a = 0.528 ± 0.001), q-SS 
(ft = 1.2, a = 0.513 ± 0.004, (3 = 0.177 ± 0.002) in isotropic 
case J|f = 0.35, TS (ft = 4.2, /3 = 1.433 ± 0.009), SS 
(ft = 4.06, q = 0.373±0.008, /3 = 1.082±0.003) in anisotropic 
case Jif = 4.0. The dotted lines have slope value 1 which 
facilitates to discern the occurence of true TS order. 



relations. The (sublinearly) diverging peak |29| of the 
momentum distribution (nfc) = r e ifcr '(bjbj +r ) oc- 

curs at A: = 7T and the quasi-condcnsation index a can 
be extracted from the formula (n n ) oc L a (a < 1). The 
calculation of (n w ) is carried out in SSE frame with al- 
gorithm of measuring Green's function 37 1. 

In Fig. [3] (a-c), we show the density peak (n v ), the 
structure factor S(tt) and the index a versus ft. In there, 
the transition points are determined by the occurring of 
meaningful nonzero a, if the fitted a < 0, that means 
there is other order instead of the triplon supcrfiuid, thus 
we denote it as zero. To check the validity of the use 
of quasi-condensation index a, we work out the transi- 
tion points at the isotropic case. For J? = J,| = 0.35, 
SSE simulation for the original ladder system gives the 
transition points from SD to SF, h c i =0.65 ± 0.03), and 
from SF to FP, h c2 = 0.88 ± 0.02) (see Fig. (c)). 
The exact solutions of it are ft c i = J± — Jii = 0.65 and 
h C 2 = J± + 2 J|| = 0.85, in good agreement with QMC. 

The triplon supersolid. — When we raise the anisotropy 
along the rail [39|, say J,f =/= J,f, the SF spreads to a 
considerable area, which can be charted by the quasi- 
condcnsation index a. If the anisotropy is large enough, 
two types of Ising orders are energetically favored. The 
first one is the traditional Ising order, i.e., AF arrange- 
ment. It is the usual Ising order and can be represented 
by IV'af) = n.j( u l s )j + ( — 1 ) J -u; | to ) j ) with u and w be- 
ing the variational parameters to minimize the MF en- 
ergy. The order parameter is the staggered magnetiza- 
tion m s = Sjp))- The transition from 
SD to AF is of second order, which is well signified by the 
divergent peak of the staggered susceptibility Xs 
(see Fig. \5ld) and (e)). 



dm s 
dh 



The second type of solid order is the alternative ar- 
rangement of singlet state |s) and triplet state \t+) on 
the rung along the rail of the ladder. It is the TS order 
represented by |^ T s) = Y\j( u A\s)2j + v A \t + ) 2j ) ® 
(ub\s)2j+i + (—I) 1 VB\t+)2j+i) with variational parame- 
ters ua ^ ug and va ^ Vb- The TS order breaks the 
translational symmetry of the ladder along the rail. To 
detect the TS order, we use the triplon structure fac- 
tor S(ir) = iE,; i (-l) ( '' +i) (^i + '»j,t + )' True TS or- 
der means a finite diagonal order parameter squared, 
("Its) 2 = r , 5'( 7r )- The true TS order is a consequence of 
the interplay between the anisotropy and the strong field. 
While in the SF regime, there may exist a weak quasi- 
triplon solid order that arises from the power-law decay- 
ing correlators j36|. It could be further strengthened by 
the anisotropy. So we define the diagonal order index /3 
that shall be extracted by the formula S(tt) oc L,p and 
set up the conditions for judging diagonal TS orders: if 
/? > 1, there exists a true TS order; if (3 < 1, a quasi-TS 
order; while if no suitable /3, no any TS order. Mean- 
while, if superfluidity (i.e., persistent LL) also survives 
under the anisotropy, then one would possibly get true 
SS order as well as q-SS order. 

We demonstrate the density peak (n w ) and the struc- 
ture factor S(tt) versus lattice size in log-log scale in Fig. 
IH From the slope of the data lines, one can figure out the 
indices a and f3. It is definite to see the coexistence of TS 
order and quasi-condensation in the SS phase. We also 
noted that the q-SS phase can occur at the isotropic case, 
as is shown in Fig. [3](c). We show the density peak (n^), 
the structure factor peak S(ir) in different lattice sizes, 
the indices a and (3 versus h (see Fig. [5]). We see the 
SS phase is a intermediate stage of the transition from 
TS to AF phase. In TS phase, the densities n s and nt+ 
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FIG. 5: (Color online) (a) The density peak (n T ) and (b) 
the structure factor S(n) versus external field h for different 
system size L. (c) The index a and j3 fitted from n{ir) and 
S(n) where the model parameters are J x = 0.35, J z = 4.0. 



prevail. Before entering into the AF phase, we find the 
system gains a considerable small fraction of n t Q <~ 0.13, 
which plays the role of impurities in the SS phase [26j ]. 
This may be the reason why SS can not emerge in general 
ID hard-core boson model or XXZ model but in this sys- 
tem. It is worthwhile to note that the simple mean-field 
description for the superfluidity in ID is not quite right. 
One can combine all the components in a wave function 
as |V) = UjNs)j + (-l) j v\t+)j + (-l) j w\t )j). But this 
assumption would lead to a true condensation and, thus, 
totally failed to capture the quasi-condensation nature. 
So the situation here is quite similar to the ID spin 1 
chain with uniaxial exchange and single-ion anisotropics, 
where the LL behavior also plays an important role in 
providing the ID superfluidity [28|. 

It is important to preclude the possibility of phase sep- 
aration. This can be done by measuring m z of the small 
area of size (say I = 8(12)) throughout the whole lat- 
tice (L = 88). If no phase separation occurs, we should 
observe only a single peak in the histograms of m z [2rj| . 
This is indeed the case as shown in Fig. El 

Experimental proposal. — The system discussed in this 
paper may be realized in a Mott insulating phase of cold 
bosonic or fcrmionic atoms in a ladder-shaped optical lat- 
tice if the atoms possess two internal states that play the 
role of the spin (such as 87 Rb in f) = \F = 1 , m p — 
1), | 1) = \F = l,mp = —1)). Let's estimate the typi- 
cal energy scales. For Rb atoms with a lattice constant 
ir/k x ~ 2ir/k y ~ 426 nm and about 10 5 atoms in a Bosc- 
Einstcin condensate, we can chose t 2 /(HU) ~ 0.1 kHz 
(corresponding to a time scale of 10 ms) with a conser- 
vative choice of U ~ 2 kHz and (t/U) 2 ~ 1/20. These 
energy scales are clearly compatible with current experi- 



ments [l[ and make the system in a Mott insulating area. 

In experiment, the density of condensates in momen- 
tum space (n„) can be measured by noise correlations 
which can be linked to spin-spin correlations [4jJ-|43j • We 
can use Bragg scattering of light, which gives rise to the 
spin structure factor, to detect S(ir) 44|. An alternative 
technique for imaging spin states in optical lattices has 
been put forward [45|. Thus, the SS discussed in this 
paper can be detected in experiment. 

Conclusions. — In summary, we propose a scheme to re- 
alize triplon SS of cold atoms in a ladder-shaped optical 
lattice. There is a q-SS phase emerging as a precursor 
of the true SS phase and TS phases. The most impor- 
tance is that this triplon SS can be detected in recent 
experimental technology These results are beyond both 
previous experimental researches in isolated double-well 
and previous studies of ladder system. It can be used as 
a guide for future experimental studies of supersolid and 
quasi-supersolid. 
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